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Abstract 

Using dispersive techniques, it is possible to avoid ultraviolet divergences in the 
calculation of Feynman diagrams, making subsequent regularization of divergent dia- 
grams unnecessary. We give a simple introduction to the most important features of 
such dispersive techniques in the framework of the so-called finite causal perturbation 
theory. The method is also applied to the 'divergent' general massive two-loop sunrise 
selfenergy diagram, where it leads directly to an analytic expression for the imaginary 
part of the diagram in accordance with the literature, whereas the real part can be 
obtained by a single integral dispersion relation. It is pointed out that dispersive 
methods have been known for decades and have been applied to several nontrivial 
Feynman diagram calculations. 
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1 Introduction 

In quantum field theory, one usually starts from fields and a Lagrangean which describes 
the interaction. These objects get quantized and S-matrix elements or Greens functions 
are then constructed by the help of Feynman rules. However, these Feynman rules lead to 
products of distributions with singular behavior. This leads to the well-known ultraviolet- 
divergences in perturbation theory. The occurrence of such divergences has led to the 
invention of several ingenious formal procedures like Pauli-Villars regularization or di- 
mensional regularization The situation remained not entirely satisfactory for several 
decades. So Feynman in his Nobel lecture remarked: "I think that the renormalization 
theory is simply a way to sweep the difficulties of the divergences of quantum electrody- 
namics under the rug." 

Today, the situation is not so bad as Feynman described it. The regularization pro- 
cedures which remove infinities from the theory are well understood and it has become 
clear that UV divergences do not imply an inconsistency of the theory in general, but 
are a consequence of our current perturbative formulation of quantum field theories. 
Renormalization- group methods also shed a new light on their interpretation. 

It is still widely unknown that it is possible to treat Feynman diagrams in such a 
manner that ultraviolet (UV) divergences do not appear in the calculations. This can 
be achieved by making explicit use of the causal structure of quantum field theory, i.e. 
by using dispersive techniques. Contrary to widespread belief, the use of dispersion rela- 
tions in practical calculations is not unusual and has been employed by various authors in 
phenomenologically important evaluations of calculationally demanding quantum electro- 
dynamic corrections such as higher-order binding corrections to the two-loop bound-state 
quantum electrodynamic self-energy. Artru et al. have used dispersive techniques as early 
as in 1966 for a leading-order calculation of the Lamb shift [SI Pachucki has employed 
these relations for the evaluation of higher-order binding corrections to the two-loop Lamb 
shift 5 , and mentioned in his paper that the ultraviolet divergences canceled automati- 
cally in the calculations. 

A very natural and rigorous approach to axiomatic perturbation theory which has 
been applied successfully to all relevant interactions of the Standard model El E] 
and supersymmetric models jl()l 111] was provided by a classic paper of H. Epstein and V. 
Glaser in 1971 |12l I13j . Their method, called finite causal perturbation theory or FCPT 
for short in this paper, avoids UV divergences from the start by defining mathematically 
correct time-ordered products for distributions. The resulting distributions are smeared 
out by test functions, avoiding provisionally volume and infrared divergences at once. An 
extensive introduction to the causal method can be found in the textbook of G. Scharf 
on finite quantum electrodynamics ^3]. There is now also a growing interest in FCPT, 
because it allows to formulate a consistent renormalization theory on a curved physical 
background [T5] . 

Basically, FCPT provides a clearly defined strategy for the dispersive calculation of 
Feynman diagrams. It renders all diagrams finite, but does not automatically solve prob- 
lems related to the specification of normalization conditions and possible anomalies arising 
in the regularization procedure, which may violate symmetries which are necessary to give 
physical meaning. These problems like e.g. the formulation of Ward identities or the 
introduction of interacting fields in the framework of FCPT have been trated extensively 
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in recent years QUI H3 U3 U3 HOI ■ 

In FCPT, the ansatz for the 5-matrix as a power series in the coupling constant is 
crucial, namely 5 is considered as a sum of smeared operator-valued distributions 

00 1 f 

S(g) = 1 + — d 4 xi . . .d 4 x n T n (xi,. . . ,x n ) g(xi) ■ . . . ■ g(x n ) , (1) 

n=l 

where the Schwartz test function g E 5(R 4 ) switches the interaction and provides the 
infrared cutoff. The basic formulation of causality in FCPT, which had already been used 
by Bogoliubov et al. is 

S(gi+g2) = S(g 2 )S( gi ) , (2) 

if the support of g<i (x) is later than suppgi in some Lorentz frame (usually denoted by 
supp52 > supp gi). The condition allows the construction of the n-point distributions T n 
as a well-defined 'renormalized' time-ordered product expressed in terms of Wick mono- 
mials of free fields : 0(x\, . . . , x n ): 

T n (x\ , . . . , Xn) — ^ ' • (D(x\j . . . , x n ) . t n (a^i X n , . . . , X n —\ %n) • (3) 

o 

where the t® are C-number distributions. T n is constructed inductively from the first or- 
der Ti(x), which describes the interaction among the quantum fields, and from the lower 
orders Tj, j = 2, . . . , n — 1 by means of Poincare covariance and causality. The inductive 
construction of the n-point distributions T n can be considered as the main technical differ- 
ence of the theory to other approaches, since all lower orders 7\, ...T n _i must be calculated 
first in order to construct T n . But the commonly chosen form of the T n 

]T e« - x° 2 ) • ... • @(x® n l - <)Ti(x ffl ) • ... • Ti^J (4) 

IT 

where the sum runs over all n! permutations, is not an unambiguous definition, since it 
contains the product of Heaviside distributions with other singular distributions. This 
leads to the well-known UV divergences in the calculation of Feynman diagrams. One 
can illustrate this fact by the following simple one-dimensional example. The product 
of a Heaviside distribution Q(t) with a ^-distribution S(t) or even its derivative S'(t) is 
obviously not well-defined. The Fourier transforms of the distributions are 

V ; ^w-iO V2T V ' 

The ill-defined product (05') (t) goes over into a non-existing convolution by Fourier trans- 
form 

{@ 5')(t) - J-(6 * 5')(u;) = / + °° dJ (6) 

which is definitely 'UV divergent'. 

It is the aim of this paper to demonstrate how UV finite results are obtained in the 
framework of FCPT by avoiding problematic Feynman integrals. This is done in Section 
2, where technical details are explained. In Section 3 the method is applied to the well- 
known introductory case of a scalar selfenergy diagram, and in section 4, we apply the 
method to the less trivial case of the two- loop sunrise diagram with arbitrary masses. 
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2 The method 

We start our considerations with a scalar neutral massive field 4> m with mass m which 
satisfies the free Klein-Gordon equation. It can be decomposed into a negative and positive 
frequency part according to 

4> m (x) = <j>m(x) + 4>ti{x) = 
(2vr)- 3 / 2 / 4i \a(k)e- ikx + a!(k)e ikx ], (7) 



'2E 

where E = \jk 2 + m? . The commutation relations for such a field are given by the Jordan- 
Pauli distributions 

[4>m( x )> <t>m{y)\ = ~iD±{x - y) , (8) 
which have the Fourier transforms 

Dtip) = (2n)- 2 J d*xD±(x)e^ = ±^G(±p°)6(p 2 - m 2 ). (9) 

Since the commutator 

[4> m (x),(j) rn (y)} = -iD~(x -y) - iD+(x - y) =: -iD m (x - y) (10) 

vanishes for spacelike distances (x — y) 2 < due to the requirement of microcausality, 
it follows immediately that the Jordan-Pauli distribution D m has causal support, i.e. it 
vanishes outside the closed forward and backward lightcone 

suppAn(x) cF~uF + , (11) 

V ± = {x\x 2 > 0, ±x° > 0}. (12) 

A crucial observation is the fact that the retarded propagator is given in real space by the 
covariant formula 

DZ\x) = 0(x°)D m (x), (13) 
which goes over into a convolution in momentum space 

t)Z\p) = ( 2 <T 2 / d 4 k b m (k)Q(p - k). (14) 

The Heaviside distribution 0(x°) could be replaced by Q(vx) with an arbitrary vector 
v £ V + inside the forward lightcone. The Fourier transform of the Heaviside distribution 
0(x°) can be easily calculated 

6(p) = (2tt)~ 2 lim / d^xQ^e-^e^ -^ = -^^6 (3) {p). (15) 

e^O J p v + z0 

For the special case where p is in the forward lightcone V + , we can go to a Lorentz frame 
where p = (p°, 0) (= p° = po as a shorthand) such that eq. (|14l) becomes 
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For arbitrary p £ V + , D™* is therefore given by the dispersion relation 

It is trivial but instructive to perform the actual calculation of Z)™*. Exploiting the 
(^-distribution in 

£> m (p) = ^-sgn(p°)5{p 2 - m 2 ), (18) 



we obtain 



1 f sgn(tp°)5(t 2 p 2 — m 2 ^ 



(2vr) 2 J 1-t + iO 



dt 



m_ 

|2 ' 



(2vr) 2 J 2y^m(l-t + i0) 

1 1 



(19) 



(27r) 2 p 2 — m 2 + iO 

We recover the analytic expression for the Feynman propagator, which coincides with 

= *Sw + dm = -j±- 2 p2 _ J - ipo0 (20) 

for p G ^ + . The imaginary part of the Feynman propagator is given by 

lm(D^p)) = ^-8(p 2 -m 2 ) (21) 

and can be deduced directly from the causal Jordan-Pauli distribution D m in an obvious 
way. 

2.1 The scalar one-loop selfenergy diagram 

As a first step, we define the causal distribution d(x — y) in conformity with (JSJ) as the 
following vacuum expectation value of a massive and a massless field (the colons denote 
normal ordering) 

H) d(x - y) := <0| [: </> m (x)</> (x) : > : <t>m(y)Mv) =]|0) (22) 
which has again causal support due to microcausality. It is sufficient to consider 
(-i) r(x - y) := -(0| : m (y)^ o (y) : : <f) m (x)(j) (x) : |0) = 

D-(x-y)D-(x-y), (23) 

since d(x — y) = r(x — y) — r(y — x). The product of the two Jordan-Pauli distributions 
in real space goes over into a convolution in momentum space 

f(p) = j^s(D-*D-)(p), (24) 
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such that d{jp) = r{p) — f{—p). 

Since we do not calculate the time- ordered product of fields, the Feynman propagators 
are replaced by Jordan-Pauli distributions. 

The integral appearing in (|24[) is readily evaluated if one exploits the ^-distribution 
stemming from the massless field contraction in 



rip) = -j-yi J d 4 q 8(-g W)e(g° - P°)S((p - qf - m 2 ). (25) 

rip) vanishes outside the closed backward lightcone due to Lorentz invariance and the two 
O-distributions in eq. (|25|). Therefore we can go to a Lorentz frame where p = (p° < 0, 0), 
and using the abbreviation E = y/iq 2 + m 2 ) leads to 

>V) = " (2^)4 / ^«(PS + VE - m 2 )6(-E - p0) = 



-^w^-^-^i&W' (26) 

and for arbitrary p we get the intermediate result 

rip) = - 3 2 7r 3 P2 ~2 m2 e (V)e(p 2 - m 2 ). (27) 

Naive use of the dispersion relation (|17|) for the real part of the sunrise diagram would 
lead to an ultraviolet divergent expression. But it can be shown ^1] that the finite part of 
the diagram is given for p £ V + by a subtracted dispersion relation (which is also called 
splitting formula) 

oo 

tretr^ _ * f m. d ^P) 



t r6t (p) = — / dt- — \p -, (28) 

w 2vr y (t-i0) w+1 (l-t + «0) V 7 

— oo 

where to is the power counting degree of the diagram, which is in our case uo = since 
the diagram is logarithmically divergent. t ret ip) is normalized such that all derivatives of 
order < uj of t ret vanish at p = 0, i.e. in the present case 

t ret (0) = 0. (29) 

These statements are only meaningful if the derivatives of order < u of t ret ip) exist in the 
sense of ordinary functions. This is assured in most cases by the existence of a massive 
field in the theory. 

We mention here that eq. (|2^)) corresponds to the dispersion relation (113,10) used in 
the familiar textbook by Landau and Lifshitz on quantum electrodynamics |22| . 

Roughly speaking, the additional term in the denominator of (|28jl has a simple expla- 
nation. Writing (|28jl for p = (p°,0) as 

oo "0 
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it becomes obvious that the division by acts as a kind of inverse differentiation 

on the causal distribution d(x) in real space. The transformed distribution has a less 
critical scaling behavior, and can be multiplied in a well-defined way by the time-ordering 
O-distribution. After the splitting, the original part of d(x) in the forward lightcone is 
recovered by differentiating uj + 1 times along the time axis, corresponding to a multiplica- 
tion of the distribution by (p°) UJ+1 in momentum space. This observation also shows the 
connection of the method to differential renormalization J2H]. For further mathematical 
details we refer also to [2*1] . 

The splitting integral can be evaluated by elementary methods 



t ret (p) 



1 



4(2vr) 



dt 



sgn(t)G(i p — m 2 )(t p — m 2 ) 
t 3 (l - t + i0)p 2 



per 



(31) 



and leads directly to 



tip) 



4(2vr) 4 



m 



p- 



log( 



m 



p 2 — iO 



+ 1 



(32) 



Distribution theory leaves the freedom to add a constant to the result, which corresponds 
to a local term c 5(x) in real space. Such local terms have to be restricted by further 
symmetry considerations in practical cases. 

We arrived thus on a very direct way to a finite expression for the scalar one-loop 
selfenergy diagram, which corresponds up to a prefactor to the regularized expression of 
the Feynman integral 

** (33) 



(k 2 + iO)((p - k) 2 - m 2 + iO) ' 

We finally mention that one-loop diagrams in two space-time dimensions are treated 
in the framework of FCPT in [251, whereas three-dimensional QED is treated in |14j . 



3 The Sunrise Diagram 

In the following section, we will utilize the strategy discussed above for the calculation of 
the imaginary part of the sunrise diagram. This will show that the causal treatment of 
Feynman diagrams is also applicable to non-trivial cases. An extensive introduction to the 
calculation of two-loop diagrams in the causal method has been given in but the case 
of the sunrise diagram is missing there. Due to the fundamental structure of the sunrise 
diagram, it is clear that it has been investigated by the use of many different approaches, 
and in this paper it serves only as a convenient example for a dispersive calculation. We 
give therefore a short overview over the literature which is related to the subject in the 
following. 

There is now an increased interest in precise calculations of multi-loop Feynman dia- 
grams. In the general massive case, relevant for future high precision calculations in the 
electroweak theory which go beyond the currently available one-loop electroweak radia- 
tive corrections |27j . the number of parameters makes it impossible to obtain results in the 
usual analytic form already in the case of the two-loop sunrise selfenergy diagram shown in 
Fig.l. When one or two internal particles are massless, the four-dimensional results can be 
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obtained in terms of dilogarithms |281 1291 I3U1 I31| . The situation is more complicated when 
all the virtual particles in the diagram have (different) masses. Such a situation occurs 
e.g. in the two-loop off-shell contributions to the Higgs selfenergy. In three dimensions, 
the expression for the sunrise diagram is quite compact even for different masses |32| . but 
in the four dimensions, there are arguments |33| that the result can not be expressed in 
terms of polylogarithms or other well-known special functions. 

There is nowadays also a widely accepted procedure of expressing radiative correction 
calculations in terms of a limited number of master integrals (MI) [H3] , which reduces the 
problem of finding analytic expressions for two-loop diagrams to the careful determination 
of the MI. The method has also the advantage that, with a correct bookkeeping of the 
recurrence relations arising from integration by parts identities, the MI of a given problem 
can be reused in more complicated calculations. 

As already mentioned, the analytical calculation of MI (in terms of the usual polylog- 
arithms and their generalizations) is possible only in cases with high symmetry, i.e. when 
the number of different scales (internal masses and external momenta or Mandelstam vari- 
ables) is small like in QCD calculations, where all masses are set to zero, or in QED-type 
cases, where only the electron mass is different from zero, or when the external variables 
are fixed to particular values (zero or mass shell condition). Another possibility of help 
in analytic calculations is sometimes offered by the exploitation of particular simplifying 
conditions, like the smallness of some ratios of the parameters allowing the corresponding 
expansion. 



m 1 




Figure I: The general massive two- loop sunrise selfenergy diagram. 

The MI of the sunrise diagram were recognized to be expressible in closed form as a 
combination of four Lauricella functions, a special class of generalized hypergeometric se- 
ries |35| (and references therein). The method provides efficient multiple series expansions 
for the regions of small \p 2 \, i.e. \p 2 \ < max(m-), and of large \p 2 \ 3> (mi + m 2 + m 3 ) 2 , 
but some problems arise in the intermediate region. 

There were also efforts devoted to the investigation of properties in special points (i.e. 
at p 2 = 0,oo, pseudothresholds and threshold). The analytical expansions of the MI 
at and oo are given in [HE] and [HE]; the values at pseudothresholds and threshold in 
[HZ]; the analytical expansions at pseudothresholds can be found in [HE]; a semi-analytical 
expansion at threshold in |39j and also in configuration space technique in [30]) while the 
complete analytical expansions at threshold are presented in |41| . 

For numerical evaluation purposes, it is possible to cast the general massive selfen- 
ergy diagram in a double integral representation and in the particular case of the sunrise 
diagram in a single integral representation |351 1421 1261 143j (and references therein) . The 
configuration space technique is also exploited in the numerical approach (HH] , [11] • In a re- 
cent approach rearrangements of the integrand, driven by the Bernstein-Tkachov theorem, 
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are introduced to improve numerical convergence JI5 . A different and interesting method 
is the use of the recurrence relations as difference equations to numerically evaluate the 

mi gnj. 

It is of general interest to provide the widest possible range of insight into the structure 
and analytic properties of multi-loop diagrams. In the present paper we exploit the causal 
structure of the sunrise diagram, which immediately leads to a straightforward evaluation 
of the imaginary part of the diagram by elliptic functions. The real part of the diagram can 
then be obtained from a single integral dispersion relation. The dispersive method works 
especially well for diagrams which depend only on one external momentum, although its 
has already been used for different kinds of cases. 

As in the case of the one-loop selfenergy diagram, we define the causal distribution 
d{x — y) according to eq. (JSJ) as the vacuum expectation value of three massive fields (the 
colons denote normal ordering) 

(-i)d(x-y) : = 
(0|[: cp mi (x)(j) m2 (x)(j) 
which has causal support. As for the one-loop case, we consider 

(-i)r(x - y) : = 

-(0| : 4>m 1 (y)4> m2 (y)4> m[i (y) ■■■ 4>m 1 (x)4> m2 (x)4> m3 (x) : |0) 

= iD mi (x - y) D m 2 ( x - y) D m 3 ( x ~ I/) > ( 35 ) 

since d(x — y) = r(x — y) — r{y — x). The product of the three Jordan-Pauli distributions 
in real space goes over into a threefold convolution in momentum space 

Hp) = -Tj^AlC^mi * D rn 2 * £>m 3 )(P)> ( 36 ) 

such that d(p) = r(p) — f(— p). 

The dispersion relation (|17|) for the real part of the sunrise diagram is given by the 
splitting formula (p £ V + ) 

iret { ) = ±_ f dt (37) 

since the power counting degree of the quadratically divergent diagram is in this case 
u) = 2. t ret (p) is normalized such that all derivatives of order < uj of t ret vanish at p = 

{".(„) -0, ^'(0> = 0, ^"*(0)=0. (38) 
3.1 The Imaginary Part 

We will now calculate the three-fold convolution of the Jordan-Pauli distributions in mo- 
mentum space given by eq. (|36|) . First we evaluate 

du(p) = (D mi *D m2 ){p) = 
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(2*) 5 



A <3>{-q°)H<l 2 - m?)e(g° - p°)<5((p - g) 2 - ml). (39) 



^12 (p) vanishes outside the closed backward lightcone due to Lorentz invariance and the 
two 0-distributions in eq. ()39|) . In a Lorentz frame where p = (p° < 0,0) we obtain 
(E=^{q 2 +m 2 )) 

di2(p°) = 

-j^yj ^(P 2 + 2 P °E + m\- m 2 )Q(-E - p°) = 
1 ( ,| i | i _ f /p° _ m 2 — 7tt-2 s 



/ d\ms(^ + E + ^^)e(-E-p°) 



2(2vr)|p°| J V 2 2p° 

-27^M 0( - /)0( ^- (mi + m2)2) 



2 2p° 

and for arbitrary p we get the intermediate result 



p^rn^^-rnl (40) 



= -^Ke(-P°MP 2 - {mi + ma) 



p 4 — 2(m 2 + m 2 ,)^ 2 + (m 2 — m 2 ,) 2 , (41) 

which is symmetric in mi, 777,2 and exhibits the typical two-particle threshold behavior. 
Applying (|28|) with u = to di2(p) would generate the retarded one-loop propagator for 
a diagram with masses mi and 777,2- As a next step we calculate 

0\r/„2 ^,2\a t (ja „0\ Q //„ „\2 / , „, \2\ 



©(-9 )% - rnl)®{q" - p°)9((p - g) 2 - (mi + m 2 

\J{p — q) A — 2(m 2 + m 2 )(p — g) 2 + (m 2 — m 2 .) 2 
(p - <?) 2 



(42) 



For p = (p°, 0) we obtain (.E = \J q 2 + 

r{p) - 



2 

m 3 



d 3 



9 



4(2vr) 6 J 2E 
Q(-E - p°)@(po + 2p E + ml - (mi + m 2 ) 2 ) 



+ 2p E + m§ 
r e(-p°)6(p§ - (mi + m 2 + ms) 2 ) 



Hpo) 



4(2vr) ; 



p | (m 1 +m 2 ) -m 3 
" 2 H 03 



/g2 _2 



/ js i?S' w ' (43) 



Pq + 2P0-E + ru- 
ms 
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I(Po) = {(p 2 + 2 Po E + m 2 ) 2 - 
2{m\ + m 2 2 )(p 2 + 2p E + m|) + (m\ - m 2 2 ) 2 Y /2 . (44) 



The integral ()43l) can be written in a compact manner if we substitute s = p 2 ) + 2poE -\-m\. 
Going back to arbitrary Lorentz frames by replacing p° by —yp^ where necessary, we 
obtain the following representation of f{p): 

Hp) = - 16(2 ^ 5p2 Q(-P°)Q(p 2 - (mi +m 2 + m 3 ) 2 ) J(p), (45) 

where 

J(P) = / -J \/( S " *l)( S " S 2)(S 3 " S)(S4 " S), (46) 

si 

and the variables s 2 < s\ < S4 < S3 are defined via 

s 1 = (m 1+ m 2 ) 2 , s 2 = ( mi -m 2 ) 2 , 

S3 = (\/p^ + m 3 ) 2 , s i = (^-m 3 ) 2 . (47) 
The integral in (|46|) can be expressed by complete elliptic integrals E, K and II: 

J = lx E(a) + l2 K{a) + 73 II(/3i, a) + 7^, a), (48) 

where 



a 



(S4 ~ glKfg - S 2 ) 
(«3 - Sl)(s4 - S 2 )' 



S4-SI a S 2 (S4 — Si) 
Pi = , P2 



S4-S2' Sl(s4-S 2 ) 



71 4 



72 = 7W — — —(si - s 2 )(si - s 2 + 3s 3 + s 4 ), 
4 V s 3 - si 

1 ( s i - s 2) fsr^ 2 n ^ 
^--4 ^ 4 -. s ,)( S3 - Sl) g'-- 2 §^»- 

74 _ _■,«*.(«.-*) (49) 

V(S4 - S2)(S3 - Si) 

For the sake of convenience and clarity, we note the definition of the complete elliptic 
integrals: 



1 



E(k) = J dt 




1-kH 2 
1-t 2 
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TT(„ k) = [ — „ dt - (50) 

1 ' ^ J ^/(l-t 2 )(l-k 2 f 2 )(l-u 2 t 2 ) V ' 

E, K and II satisfy very many identities, therefore it might be possible that an even more 
compact form can be found for J. The analytic expression for d(p) follows from f(p) simply 
by replacing 0(— p°) by — sgn(j>°), and the imaginary part of the sunrise diagram is given 
by 

Im = 32(2 ^ 5p2 ©CP 2 ~ (mi + m 2 + m 3 ) 2 ) J . (51) 



The result is indeed fully symmetric under permutations of mi, 777-2,7113, although this is 
not obvious from eq. (|48|) . But for the higher symmetry cases, J becomes quite simple. 
For mi = 0, 7772 = m 3 = 777, J reduces to 



J(p) = I(p 2 + 2m 2 )^p 2 (p 2 - 4m 2 ), (52) 



for mi = w-2 = 0, 7773 = 777 we have 

1 



J(P) = ^(j7 2 - 777 2 )(j7 2 +m 2 ), (53) 

and for mi = m 2 = 7773 = , J collapses to 

J(p) = i/. (54) 

Applying eq. (|3*7j) to (|52*|) and (fB^|) leads to the results mentioned in the introduction which 
are expressible by polylogarithms and associated functions. But the subtracted dispersion 
relation (|37|) does not work in the totally massless case, where one encounters an additional 
infrared divergence which requires additional regularization that can be derived also as a 
limit from the massive case. The solution in the massless case, relevant e.g. for QCD or 
quantum gravity, is 

t™'(p 2 ) = t(p 2 )+f(p) = 

1 r/log^^, (55) 



32(2tt) 6 ^ to M 2 

t( 2 \ 1 2, -(p 2 + i0) 

^ ) = "32W P bg M 2 



r^7TP 2 (log-^ -^(p 2 )) (56) 

32(2vr) 6 ^ V M 2 yF ' ) v ' 

with a real scaling parameter M. Eq. (|56|) reflects also the high energy behavior of the 
massive propagator t ~ p 2 log(p 2 ). The prefactor (2?t)~ 6 in the expressions above is due 
to the fact that we have included a prefactor (27r) -2 in our definition of the Feynman 
propagator. 
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4 Final remarks and conclusion 

A closed analytic expression has been derived for the imaginary part of the scalar sunrise 
diagram with arbitrary masses of the virtual particles. The technical difference between 
the usual Feynman diagram calculations and the method used in the present paper is 
that the occurring integrals are reduced step by step to single integrals by utilizing the 
^-distributions in the Jordan-Pauli distributions instead of performing a Wick rotation. 
The real part of the sunrise diagram can be expressed for all values of the external impulse 
by a single integral, because the causal distribution d which is related to the imaginary 
part of the diagram has the special form 

d{p) = /i(p 2 )sgn(> )eO 2 - mf ot ), mtot = m 1 + m 2 + m 3 , (57) 

therefore the dispersion relation (|37|) can also be written for p £ V + 

f ret (p) = t(p) + rip) = 

oo 

m tot 

where [uj/2] denotes the largest integer <u/2. This formula can be extended to space-like 
p by analytic continuation 

oo 

m tot 

such that that an access to numerical and analytic investigations is provided for the whole 
range of external momenta. 

The problem of overlapping divergences is automatically solved by the causal method 
due to the inductive construction of the theory, whereas in usual renormalization schemes 
it becomes a nontrivial problem to show that all sub-diagrams of a given diagram can be 
renormalized in a consistent way. In the context of the BPHZ method, the problem was 
solved by the forest formula, an algorithm which disentangles all divergences in subdia- 
grams This is one of the conceptual strengths of the causal approach. Furthermore it 
should be pointed out that on the level of the sunrise diagram discussed in this paper, the 
calculation of the imaginary part can as well be understood by the Cutkoski rules. But 
at higher orders, one should clearly distinguish between these rules which rely mainly on 
the unitarity of the theory, whereas in the causal approach the main input is causality. 

We have used the example of the surise diagram in order to illustrate in a pedagogical 
manner the strategy and fundamental properties of the dispersive approach for the cal- 
culation of Feynman diagrams. We refer also to the recent literature on the specific case 
of the sunrise diagram, where also non-trivial vertex structures are treated by dispersive 
and non-dispersive methods [J2l HU EH] • 
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